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Abstract 

Background: Maternal undernutrition leads to an increased risk of metabolic disorders in offspring including 
obesity and insulin resistance, thought to be due to a programmed thrifty phenotype which is inappropriate for a 
subsequent richer nutritional environment. In a rat model, both male and female offspring of undernourished 
mothers are programmed to become obese, however postnatal leptin treatment gives discordant results between 
males and females. Leptin treatment is able to rescue the adverse programming effects in the female offspring of 
undernourished mothers, but not in their male offspring. Additionally, in these rats, postnatal leptin treatment of 
offspring from normally-nourished mothers programmes their male offspring to develop obesity in later life, while 
there is no comparable effect in their female offspring. 

Results: We show by microarray analysis of the female liver transcriptome that both maternal undernutrition and 
postnatal leptin treatment independently induce a similar thrifty transcriptional programme affecting carbohydrate 
metabolism, amino acid metabolism and oxidative stress genes. Paradoxically, however, the combination of both 
stimuli restores a more normal transcriptional environment. This demonstrates that "leptin reversal" is a global 
phenomenon affecting all genes involved in fetal programming by maternal undernourishment and leptin treatment. 
The thrifty transcriptional programme was associated with pro-inflammatory markers and downregulation of adaptive 
immune mediators, particularly MHC class I genes, suggesting a deficit in antigen presentation in these offspring. 

Conclusions: We propose a revised model of developmental programming reconciling the male and female 
observations, in which there are two competing programmes which collectively drive liver transcription. The first 
element is a thrifty metabolic phenotype induced by early life growth restriction independently of leptin levels. The 
second is a homeostatic set point calibrated in response to postnatal leptin surge, which is able to over-ride the 
metabolic programme. This "calibration model" for the postnatal leptin surge, if applicable in humans, may have 
implications for understanding responses to catch-up growth in infants. Additionally, the identification of an antigen 
presentation deficit associated with metabolic thriftiness may relate to a previously observed correlation between birth 
season (a proxy for gestational undernutrition) and infectious disease mortality in rural African communities. 
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Background 

Alterations in nutrition during fetal and perinatal life are 
linked to adverse health outcomes in offspring in adult- 
hood, this being known as the Developmental Origins of 
Health and Disease (DOHaD) paradigm [1]. In particu- 
lar, data from epidemiological cohorts and a range of 
animal models has shown that maternal undernutrition 
results in an increased risk of obesity and metabolic dis- 
ease in offspring in later life [2,3]. The most commonly 
used model of maternal undernutrition is in the rat where 
studies have utilised moderate through to severe undernu- 
trition (20-70% calorie restricted diets) to examine mecha- 
nisms underlying the programming of later disease risk 
[4-6]. Protein undernutrition appears to be particularly 
important, with low protein isocaloric diet models show- 
ing similar effects to caloric restriction models in both rat 
and mouse [7,8]. 

The adverse phenotypic outcomes, particularly those 
related to metabolic abnormalities, are thought to result 
in part from the mismatch between a deprived early nutri- 
tional environment which programmes a "thrifty phenotype" 
(also called a "predictive adaptive response"), and a richer 
later environment for which this metabolic programme is 
inappropriate [9-11]. In addition to the metabolic sequelae 
of maternal undernutrition, in some circumstances such 
deprivation can have long-lasting consequences for the 
immune system, for example the greatly increased preva- 
lence of infectious disease among Gambian individuals 
born in the "hungry season" [12-14]. This latter result was 
however not replicated in studies in rural Bangladesh or 
Senegal [15,16], indicating a wide degree of heterogeneity 
in the immune response to maternal undernutrition. 

We have previously shown in a rat model (70% caloric 
restriction) that neonatal administration of the adipokine 
leptin reverses the metabolic abnormalities seen in the 
offspring of undernourished mothers [17,18]. Both male 
and female offspring of mothers subject to maternal under- 
nutrition during pregnancy developed increased adiposity 
and markers of the metabolic syndrome, particularly when 
fed a high fat diet postweaning. Neonatal treatment with 
leptin during a key period of developmental plasticity 
(postnatal age 3-13 days) reversed the programmed pheno- 
type and restored near normal metabolic parameters 
[17,18]. This observation was associated with a rever- 
sal of the direction of the leptin response (hereafter "leptin 
reversal") in the offspring of undernourished mothers 
('UN offspring'). Neonatal leptin treatment increased 
adult levels of llp-HSD2 transcription in the livers of 
normally-nourished pups, but decreased it in livers of the 
UN offspring: similar leptin reversal effects on both tran- 
scription and promoter methylation were observed for 
PPARa, GR and PEPCK [19]. 

In this follow-up analysis, we carried out whole- 
genome expression profiling of liver RNA from the female 



experimental series previously described, in order to fully 
characterise the global hepatic responses to maternal un- 
dernutrition (AD /UN = ad-libitum fed or undernourished 
mothers), postnatal leptin administration (Lep/Sal = leptin 
or saline control given postnatally), level of postweaning 
diet (Chow/HF = normal or high fat diet), and the interac- 
tions between these three factors. In particular, we aimed 
to determine: 

1. The extent of the leptin reversal- whether it affected 
all leptin-regulated genes, or a specific subset of 
these. 

2. Which of the genes affected by the leptin reversal 
were associated with the metabolic syndrome seen 
in UN/Sal/HF offspring. 

3. Whether there were precursor changes seen in 
UN/Sal/Chow offspring, which could potentially be 
markers for those at risk of developing metabolic 
syndrome and/or diabetes. 

4. Whether there were other non-metabolic genes and 
pathways affected by maternal undernutrition 
and/or leptin reversal which might therefore be 
implicated in further phenotypes, We were particularly 
interested in pathways relating to immune system 
regulation, given the programmed effects (in human) 
of undernutrition on infectious disease prevalence. 

The transcriptional data indicated a significant similarity 
between the AD/Lep cohorts and the UN/Sal cohorts. 
This was unexpected since the phenotypic outcomes 
differ markedly between these two interventions, in a 
sex-specific manner. In particular, in females UN/Sal 
offspring become severely obese when fed a high fat 
diet, while AD/Lep do not: the reverse is the case in 
males. These transcriptional findings therefore led us to 
re-examine our previously published growth curves for 
both the male and female experimental series. Previously, 
we had analysed the preweaning and postweaning data 
separately. In this re-analysis we focus on the weaning 
period itself: a critical juncture during which growing 
pups take over full responsibility for their own nutritional 
intake. We find that growth trajectories differ between 
male and female pups during weaning, and interpret the 
male and female data collectively as the consequences of 
two competing programmes: a thriftiness programme 
which governs the efficiency of fuel, and a homeostatic set 
point governing body composition. 

Results 

We used Illumina oligonucleotide arrays to perform ex- 
pression profiling on RNA extracted from livers of fe- 
male rats in 8 treatment groups (i.e. all combinations of 
ad-libitum-fed or undernourished mothers, postnatal lep- 
tin treatment or saline, postweaning high fat or normal 
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chow diet), n = 8 animals per group. An initial one-way 
ANOVA analysis was used to select all genes significantly 
differentially expressed in at least one of the treatment 
groups. 2221/8497 (26.1%) of liver-expressed genes were 
called as significant in this analysis, with 1069 of these 
showing at least 1.25 fold change in transcript abundance 
between the highest-and lowest-expressing treatment 
groups. Three-factor ANOVA analysis was subsequently 
used to categorise the genes according to which individual 
experimental factors, or interactions between factors, were 
significant in each case (Table 1 and Additional file 1: 
Table SI). 

Maternal undernutrition and postnatal leptin treatment 
predominantly affect the same target set 

Whether considering all significant transcripts, or only 
those showing more than 1.25 fold change, the largest 
group of transcriptional changes (category C) represents 
the non- interacting effects of postweaning diet-unsur- 
prising as diet at sacrifice directly affects the animals' 
metabolic status and consequently hepatic expression 
profiles. However, category AB (significant regulation by 
the interaction between maternal nutrition and leptin 

Table 1 Numbers of significantly regulated genes, 
categorised according to which factors and/or 
interactions were significant in the three-way ANOVA 



Category Significant Biological interpretation 

genes 



A Maternal nutrition 


368 (151) 


Maternal nutritional status has an 
effect independent of leptin 
treatment and postweaning diet 


B Leptin treatment 


200 (90) 


Leptin treatment has an effect 
independent of maternal nutritional 
status and postweaning diet 


C Postweaning diet 


990 (531) 


Postweaning diet has an effect 
independent of maternal nutritional 
status and leptin treatment 


AB (2-way interaction) 


949 (503) 


The effect of leptin treatment is 
dependent on maternal nutritional 
status (or vice versa) 


AC (2-way interaction) 


96 (18) 


The effect of postweaning diet is 
dependent on maternal nutritional 
status (or vice versa) 


BC (2-way interaction) 


127 (74) 


The effect of postweaning diet is 
dependent on leptin treatment 
(or vice versa) 


ABC (3-way interaction) 


7(3) 


The effects of all three factors are 
mutually dependent 


Uncategorised 


128 (59) 


Although significant changes were 
detected by one-way ANOVA, no 
individual factor or interaction term 
was subsequently significant in a 
three-factor ANOVA 



Bracketed figures for each category indicate the number of genes in the 
category that showed at least 1.25 fold change in transcript abundance across 
the data set as a whole. Additional file 2: Table S2 shows expression data for 
all genes falling into each category. 



treatment) contained very nearly the same number of 
transcripts, indicating that (a) these two factors are able to 
collectively programme transcriptional responses occur- 
ring many days post-intervention, and (b) these program- 
ming effects are equally as important as postweaning diet 
in determining the final liver expression profile. Import- 
antly, category AB contains many more genes than either 
category A or B alone, indicating that maternal undernu- 
trition and postnatal leptin treatment largely affect the 
same genes and that these effects are mutually dependent. 

Phenotypic data for the programmed cohorts also show 
significant interactions between maternal diet and 
postnatal leptin treatment 

The phenotypic characterisation of this experimental 
series was expanded and re-examined, again using three- 
factor ANOVA (Additional file 2: Table S2). Importantly, 
while total body fat percentage is increased as expected 
by a postweaning high fat diet (postweaning diet F 1)56 = 
150.5, p = 1.67 x 10" 17 ), it also shows a significant inter- 
action between maternal undernutrition and postnatal 
leptin treatment (AB interaction F 1>56 = 22.25, p = 1.63 x 
10" 5 ). Total body fat percentage is increased by postnatal 
leptin treatment (AD/Lep) and maternal undernutrition 
(UN/Sal) but decreased by the combination of both (UN/ 
Lep). Plasma leptin levels follow the same pattern as body 
fat percent (AB interaction F 1>56 = 11.32, p = 1.39 x 10" 3 ), 
and fasting C-peptide levels also show a highly signifi- 
cant AB interaction term (AB interaction F 1>56 = 13.13, 
p = 6.29 x 10" 4 ) and a similar profile of dysregulation. 
The AB interaction term for fasting insulin levels was just 
significant before FDR correction but not significant after 
correction. Plasma ghrelin and other metabolic parame- 
ters (liver glycerol and triglycerides, plasma glycerol, tri- 
glycerides and free fatty acids) were all predominantly 
regulated by postweaning diet. 

There is a complex spectrum of "reversal" interactions 
between transcriptional programming by maternal 
nutritional status and by postnatal leptin treatment 

The interactions between maternal nutrition status and 
leptin treatment could be either synergistic or opposing. 
For the 949 genes in category AB, we compared the 
expression change induced by leptin treatment in pups 
born to UN mothers (average Lep/Sal expression ratio av- 
eraged across UN cohorts) to that in pups of ad libitum 
fed mothers (the same ratio averaged across AD cohorts) 
to determine the interaction type. The large majority of 
genes (871/949 = 91.8%) showed a "leptin reversal" pattern 
where the average Lep/Sal ratio in UN cohorts was oppos- 
ite in sign to that in AD cohorts. Only a small proportion 
of genes showed a synergistic interaction (57/949 = 6.0%) 
or a partially opposing interaction (21/949 = 2.2%), defined 
as genes where the average Lep/Sal ratio in UN cohorts 
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was the same sign as in AD cohorts and of higher magni- 
tude (synergistic) or lower magnitude (partially opposing). 

Hierarchical clustering of the category AB genes was 
performed to determine whether there was evidence for 
a binary polyphenism as suggested by Gluckman et al 
[19]. Figure 1 is a heatmap of these 949 genes showing the 
residual expression changes attributable to the combin- 
ation of maternal diet and/or leptin treatment. Interacting 
genes fell into six broad groups of response pattern, all of 
which show "leptin reversal" since the direction of leptin- 
related change is opposite in AD and UN cohorts, but 
with important differences in the nature of the reversal. 

Group 1 and 3 genes are induced both by maternal 
undernutrition and by leptin treatment (compare USC/ 
H and ALC/H to ASC/H), however the combination of 
both factors (compare ULC/H to ASC/H) causes either 
no induction of these genes (Group 3) or even a net re- 
pression (Group 1). Group 2 genes show no net effect of 
maternal undernutrition on its own (compare USC/H to 
ASC/H), but are induced by leptin in AD cohorts (com- 
pare ALC/H to ASC/H) and repressed by leptin in UN 
cohorts (compare ULC/H to USC/H). Group 4 and 5 
genes are repressed by both factors acting alone, but 
show a reduced repression or even a net induction when 
both factors are present (opposite to groups 1 and 3). 
Group 6 genes are induced specifically by the combin- 
ation of both factors. 



The vast majority of AB-interacting genes show similar 
changes in AD/Lep cohorts and in UN/Sal cohorts rela- 
tive to the control AD/Sal cohorts. These results there- 
fore support the hypothesis of a binary polyphenism in 
liver transcriptional activity, where the phenotypic switch 
is seen in response both to undernutrition and to postna- 
tal leptin treatment, but not the combination of both. 

Functional annotation analysis of leptin reversal genes 

DAVID [20] was used to collate functional annotations 
into annotation groups associated with each of the identi- 
fied AB-interacting gene clusters by combining data from 
multiple databases including Gene Ontology, Swiss- Prot, 
Uni-Prot, Protein Information Resource, InterPro and 
KEGG pathways (Table 2 and Additional file 3: Table S3). 
Each annotation group has an enrichment score denoting 
whether it is significantly associated with any given gene 
cluster (enrichment score of>2 represents a p-value< 
0.01 after correction for multiple testing). This analysis re- 
vealed several functional shifts associated with maternal 
undernutrition and leptin treatment. 

(a) A series of metabolic and regulatory pathways were 
significantly over-represented in groups 1 and 3 
(i.e. upregulated in AD/Lep and UN/Sal livers), 
particularly involving mitochondrially targeted genes. 
The observed shifts are consistent with thrifty 
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Figure 1 Heatmap showing hierarchical clustering of the 949 genes with a significant interaction between maternal diet and leptin 
treatment after subtracting out the effect of postweaning diet. Line graphs show the mean (centroid) expression profile for each cluster. 
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Table 2 Functional clusters over-represented among genes with a significant interaction between maternal 
undernutrition and postnatal leptin treatment 



Expression pattern 

Induced by maternal undernutrition 
and by leptin treatment, but not 
by both combined (Groups 1 & 3) 



Annotation group 

Mitochondrially-targeted genes 
(enrichment score 2.64— » 4.55) 



Microsome/peroxisome/cytochrome 
p450 pathway (enrichment score 
2.03 -> 4.46) 



Steroid synthesis and response 
(enrichment score 1 .86 — > 2.69) 



Repressed by leptin in UN offspring, 
induced by leptin in control 
offspring (Group2) 



Starch metabolism (enrichment 
score 3.09) 

Monosaccharide metabolism 
(enrichment score 2.84) 



Pyruvate metabolism (enrichment 
score 2.23) 

Complement cascade and innate 
immunity (enrichment score 0.98) 

Metal ion binding (enrichment 
score 3.00) 



Repressed by maternal undernutrition MHC antigens and lymphocyte 



and by leptin treatment, but not 
by both combined (Groups 4 & 5) 



Induced by the combination of 
maternal undernutrition and leptin 
treatment, but not by either 
factor alone (Group 6) 



activation (enrichment score 
2.43 -> 4.00) 

Ribosomal biogenesis/translational 
elongation (enrichment score 2.27) 



Ribosomal biogenesis/nucleolus 
(enrichment score 4.87) 



Cell cycle/Mitosis (enrichment 
score 1 .70) 



Genes with AB interaction 

Aadac; Aadat; Abcbll; Acat2; Acsl3; Acsl5; Adhfel; Agtrla; Aifml; Ak3l1; 
Aldhlbl; Aldh1l2; Atadl; Atp6vlal; Bche; Bckdhb; Bphl; Cabcl; Clpx; 
Crlsl; Cry}; Cyb5r3; Cyp2r1; Dhcr24; Dhtkdl; Dnaja3; Fam82a; Fut8; 
Gludl; Gm2a; Golph3; Gpam; Gpd2; Hmbs; Hmgcr; NIT2; LOC501346; 
LOC685778; LOC688587; TTC19; Lrp5; Lyplal; Maoa; Maob; Mel; Mgstl; 
Mpv17l; Mrpl18; Mrpl50; Mterfd3; Mtfrl; Nnt; Nsdhl; Otc; Pck2; Pdhb; 
Qdpr; Rab5a; RGD13081 14; RGD1309676; PEX1; MY06; MUT; TMEM70; 
Rilp; RpllOo; Sors2; Scp2; Semi; Sec22b; Senp2; Sfxn5; Slc16o1; SI c2 5a 17; 
Slc35a3; Slc35b3; Sod2; St3gal6; Synj2bp; Tbc1d15; Tm7sf2; Tmpo; Tnks2; 
Txndd; Uros 

A lef; Aadac; Abca6; Abcb 1 I; Acsl3; Acsl5; Agl; Aifm I; Aldh 1a1; Boat; 
Bche; Cabcl; Cadml; Clcn2; Cnbp; Csnk1g3; Ctsc; Cyb5r3; Cyp2d22; 
Cyp2d5; Cyp2r1; Dhcr24; Dio1; Dpys; Ephxl; Erapl; Fmr1; Gbel; Hmgcr; 
Isod; Maoa; Mbl1; Mbl2; Mel; Mgstl; Mpv17l; Nsdhl; Nt5e; Pck2; Pgd; 
Pgk1; Pipox; Ppp1r3b; Ptprf; Pygl; Rab5a; Rdh3; PEX1; GNPNAT1; Rgs16; 
Sat1; Scp2; Sec22b; SlclOal; Slc25a17; Sod2; Sord; Srd5a1; Tm7sf2; 
Tmcol; Tmed5; Tnks2; Txndd; Ugt2b; Ugt2b17; Ugt2b36 

Aadac; Abcbll; Abhd5; Acsl3; Acsl5; Afp; Agl; Agtrla; Ak3l1; Aldhlal; 
Angptl4; Avprla; Baat; Bche; Bckdhb; Btgl; Cfb; Crlsl; Ctsc; Cyb5r3; 
Cyp2d22; Cyp2rl; Dhcr24; Dnaja3; Dnajb5; Ephxl; Erapl; Foxola; Fut8; 
Gpam; Hmgcr; Ihh; Laptm4b; LOC501346; ANGPTL3; LOC686548; 
LOC688587; Lrp5; Mel; Mgstl; Mif; Nfe2l2; Nsdhl; Osbp; Pck2; Pggtlb; 
Prkaal; Rab5a; Rcanl; RGDI5605I3; MY06; STARD5; Scrnl; Sec22b; 
Senp2; Serpina3m; Serpinfl; Slc34a2; Slc35a3; Smad2; Sord; Srd5al; 
St3gal6; Sultlal; Tm7sf2; Tmpo; Tnks2; Txndd; Uros 

Adh6; Agl; Aldhlal; Aldhlbl; Amyla; Cyp2d22; Cyp2d5; Dpys; Ephxl; 
Fut8; Gbel; Hmbs; Maoa; Maob; Mgstl; Pygl; Rdh3; Srd5al; St3gal6; 
Ugt2a3; Ugt2b; Ugt2bl0; Ugt2bl7; Ugt2b36; Uros 

Adh6; Agl; Aifml; Aldhlal; Aldhlbl; Atf4; Atp6vlal; Bckdhb; Csnklg3; 
Cyb5r3; Dhtkdl; Diol; Dpys; Gbel; Gm2a; Gpd2; Hmgcr; LOC685778; 
Maoa; Mel; Pck2; Pdhb; Pgd; Pgkl; Ppplr3b; Prkaal; Pygl; UAP1; 
GNPNAT1; Satl; SlcWal; Sod2; Sord 

Acat2; Acotl2; Adh6; Aldhlbl; Atf4; Gpd2; LOC685778; Mel; Pck2; Pdhb; 
Pgkl; Sid 6al 

Baat; C2; Cabcl; Cadml; Cfb; Dhcr24; Dnaja3; Erapl; F9; Ihh; LEAP2; 
LOC686548; Mbll; Mbl2; Mif; Mstl; Nrep; RGD1560513; Serpina5; Sod2 

Aoxl; Arsk; Atp2bl; Cdol; Cnot6; Cyp2b3; Cyp3a3; Panel; Gatad2b; Ireb2; 
ZZZ3; PHF20L1; Mdm4; Msl2ll; Nr5a2; Pcgf4; Pde8a; Phospho2; Prickle! ; 
RGD1563633; Slcl0a5; Slc30a7; Slc40al; Thrb; Tmlhe; Trim33; Upbl; 
Zfpl31;Zfp410;Zfp99 

Cc/5; Ccrl2; Cd69; Cd8a; Corola; Ctsll; Fas; Gpsm3; Gstpl; H2-M3; Icam2; 
Igh-la; lllb; Itgb2; Klrdl; RGD1561628; Msh2; P2ryl4; SH2D1A; Rpl9; RT1- 
149; RT1-A2; RT1-A3; RTl-Ba; RT1-CE15; RT1-CE7; RT1-M6-2; Sartl; Spn; 
Tap2; Tnfaip8l2; Tnfsf13; Uncl3d; Unc93bl; Vcaml; Wdr46; Xcll 

Anxal; Blocls2; Cpsf3; Epb4.1 II; Gtf2b; Hnrpab; RPL22L2; Mrps25; 
Nufipl; Nup210; Polr2h; Polr2i; Polr2j; Pop4; RGD1559639; RGD1 559951; 
Rnasen; Rpll7; Rpll8a; Rpl9; Rpsl4; Rpsl6; Rpsl7; Rps27; Sartl; Snrplc; 
Tdrd3; Tubb5; UtpHa; Uxt; Wdr46 

Agxt; Cct2; Cdc25a; Cdkl05; Colq; Copbl; Dctn6; Epb4.1 II; Fabp2; Ft/ 7; 
Gabarap; Gtf2a2; Gys2; Hnrnpal; Hspbl; Imp3; KrtlO; RPL22L1; 
LOC364236; RGD1 563484; RPS7; RGD1 '564519; RGD1 564606; POLR2E; 
Mapre3; Med6; Mettll la; MGC1 14381; Mrpl49; Myo5b; My rip; Nbn; Pola2; 
Psma6; Psma7; Psmb2; Ptk2; Ranbpl; Rassf5; WDR83; RGD1559574; 
RGD1559846; RGD1561086; RPS19L1; WTAP; RGD1565170; ZCCHC17; 
RGDI 566137; RGDI 566326; Rpain; Rpl27; Rpl35a; Rslldl; Sfrs2; Sqstml; 
Tp53; Tuba4a; Vps33b 

Cdc25a; Gadd45gipl; Mapre3; MGC1 14381; Nbn; Psma6; Psma7; Psmb2; 
Ranbpl; RgslO; Tp53; Uspl6 



Note that individual genes may be annotated within several different functional clusters. Full DAVID output for each of the six groups is given in Additional file 3: 
Table S3. 
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metabolic reprogramming in these cohorts, and 
included several key enzymes involved in glycogenolysis 
and glycolysis, such as glygogen debranching enzyme 
(Agl), glycogen phosphorylase (Pygl), phosphoglycerate 
kinase (Pgk) and two subunits of pyruvate 
dehydrogenase (Pdhb and LOC685778). Also 
upregulated were important enzymes for amino 
acid and nitrogen metabolism, including glutamate 
dehydrogenase (Glud) and ornithine transcarbamoylase 
(Otc); lipid metabolism including glycerol phosphate 
dehydrogenase (Gpd2), long-chain fatty acid-CoA 
ligases (Acsl3 and AcslS) and acetyl- Co A 
acetyltransferase (Acat2); and oxidative stress 
genes including phosphogluconate dehydrogenase 
(Pgd) and mitochondrial outer membrane 
glutathione S -transferase (Mgstl). Consistent with 
the latter, although it did not fall into any of the 
annotated functional groups in the DAVID analysis, 
carbonic anhydrase 3 (Ca3) was also upregulated in 
AD/Lep and UN/Sal livers. 

(b) Immune-related pathways showed a concerted 
pattern of change across the data set. Groups 4 and 

5 (i.e. downregulated in AD/Lep and UN/Sal livers) 
were enriched for genes involved in the adaptive 
immune response, antigen presentation and 
lymphocyte activation. Particularly striking was a 
strong downregulation of MHC genes (-4-8 fold 
change of several class I genes, but also some class 
lb and class II genes), and also a slight 
downregulation of the Tap2 peptide transporter 
necessary for antigen presentation by class I 
molecules. In this light it is interesting to note that 
groups 1-3 (with the opposite expression pattern) 
showed a trend towards enrichment for genes 
involved in the complement cascade and innate 
immunity such as complement gene C2. However, 
this functional group was not statistically significant 
(enrichment score of 0.98, corresponding to 
p-value = 0.10). We note that alpha-fetoprotein 
(Afp), a specific marker for liver inflammation also 
falls into group 3, although it was not annotated as 
inflammation-related in the DAVID analysis. 

(c) Pathways associated with ribosome biogenesis and 
function were significantly over-represented in group 

6 and (to a lesser extent) groups 4 and 5. Finally, 
mitotic cell cycle genes showed weak over-represention 
in group 6. The significance level of this was borderline, 
with an enrichment score of 1.70 (p-value = 0.02). 

Collectively, these obervations indicate that there is a 
programmed thrifty liver phenotype in AD/Lep and UN/ 
Sal livers that includes increased glycogenolysis, an in- 
crease in protein breakdown and fat synthesis, and a 
switch away from adaptive immune function towards 



innate immunity/inflammation, together with increased 
oxidative stress within the cell. 

"Healthily thrifty" versus "unhealthily thrifty": is there a 
difference? 

Our finding that both AD/Lep and UN/Sal treatments 
appear to have similar transcriptional effects was unex- 
pected, since in female rats only UN/Sal/HF develops 
metabolic syndrome while AD/Lep/HF does not. We 
therefore considered it possible that AD/Lep might rep- 
resent a state of "healthy thriftiness"-i.e. a thrifty metab- 
olism that is nevertheless under control and able to cope 
with a high fat diet without developing obesity; and that 
UN/Sal might represent "unhealthy thriftiness"-i.e. a 
thrifty metabolism that lacks safeguards and runs out of 
control when confronted with a high fat diet. 

To test this possibility, we carried out a pairwise com- 
parison of the AD/Lep/Chow and UN/Sal/Chow cohorts 
to see if there were any expression differences which 
might explain why the latter is predisposed to metabolic 
syndrome. However, we found no genes that were sig- 
nificantly differentially expressed (FDR p-value < 0.05) 
between these two cohorts, suggesting that the differ- 
ence is one of degree rather than kind. Consistent with 
this, of 61 genes significantly differentially expressed (FDR 
p-value < 0.05) between UN/Sal/Chow and the control 
AD/Sal/Chow cohort, 58 showed the same direction of 
change in AD/Lep/Chow. 

Is thriftiness a consequence of growth restriction rather 
than leptin signalling per se? 

It is challenging to explain why the UN/ Sal and AD/Lep 
cohorts both show a thrifty transcriptional profile, given 
that postnatal leptin administration in the latter case will 
mimic excess nutrition rather than undernutrition. One 
possible explanation is that the dose/response curve for 
leptin is intrinsically bathtub-shaped, with adverse ef- 
fects mediated both by excess and by insufficient leptin 
stimulation. This would be consistent with work showing 
that metabolic syndrome can be triggered by maternal 
overnutrition as well as maternal undernutrition [21]. How- 
ever, an alternative and in our view more attractive explan- 
ation is that the commonalities in expression profiles for 
these cohorts reflect a shared early life history of growth 
retardation. In the case of the UN/ Sal cohorts, the growth 
retardation is a direct consequence of the maternal under- 
nutrition, while in the case of AD/Lep cohorts it occurs 
postnatally and is a secondary consequence of the leptin 
treatment. 

Figure 2 A shows that while UN/ Sal pups are severely 
growth restricted in utero as a consequence of the ma- 
ternal undernutrition, AD/Lep pups also show a signifi- 
cant restriction in neonatal growth rate during the period 
of treatment (see Additional file 4: Table S4 for raw data). 
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Although transient, the growth deficit is sufficient to in- 
duce a 5.9 percentage point drop in body mass relative to 
control, with AD/Lep body mass being 95.7% of AD/Sal at 
day 3, and 89.8% of AD/Sal at day 14. Other studies have 
shown that growth restriction in leptin-treated neonates is 
mediated by non-hepatic effects of leptin, particularly 
increased thermogenesis [22]. In the female experimental 
series analysed here, although the leptin-treated animals 
showed no decrease in food intake relative to saline- 
treated animals during the neonatal period and are thus 
not hypophagic [17], nevertheless it indicates that they 
have a negative energy balance given their increased ener- 
getic requirements and fail to take in sufficient nutrition 
to maintain a normal growth rate. 

The leptin rescue paradox: a conflict between thriftiness 
and set point programming 

If maternal undernutrition and postnatal leptin treat- 
ment can both reprogramme liver gene expression to be 
thrifty, why then does the combination of both not do 
so? One possibility is that the postnatal leptin treatment 
might expunge the programme set during fetal life and 
then somehow fail to impose the postnatal programming 
seen in AD/Lep. Figure 2B shows that the postnatal lep- 
tin administration has the same growth- retarding effects 



in UN/Lep as in AD/Lep, rendering this hypothesis 
unlikely. In fact, the degree of growth suppression is con- 
siderably greater in the UN group, suggesting that the ef- 
fect of leptin is potentiated by the prior starvation 
episode. If the growth restriction seen in AD/Lep is suffi- 
cient to trigger hepatic reprogramming, then it is reason- 
able to assume that the even greater effect in UN/Lep 
should also leave the liver in a thrifty state. 

Figure 2B also shows that there is a significant drop 
in relative growth of the UN/Lep cohorts around wean- 
ing, which first becomes visible in the 1-2 days immedi- 
ately prior to separation from the dam (i.e. when the 
pups begin to self-wean), and continues throughout 
the following week. This juncture marks the period in 
which pups take over full responsibility for their own 
metabolic balance between intake and expenditure, ra- 
ther than it being partly controlled by maternal milk 
availability and composition. This postweaning drop 
in the UN/Lep cohort is not simply due to increased 
catch-up growth in UN/Sal, but reflects a genuine 
drop in weight relative to AD/Sal (Figure 2C). These 
observations support an explanation for the apparent 
paradox in which the two growth restriction episodes 
(both fetal and postnatal) in the UN/Lep cohorts do 
lead to thrifty metabolic programming in the livers of 
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this group, but this is masked by non-hepatic pro- 
gramming effects of leptin on food intake and energy 
expenditure. 

Re-interpreting existing male growth data in the light of 
the female transcriptional findings 

In the light of the above findings, we also re-examined 
the growth data from our previously published experi- 
ment on male rats ([18] and see also Additional file 4: 
Table S4 for raw data). Females and males show consid- 
erable differences in their programmed growth responses 
to postnatal leptin administration as early as weaning, 
where in males there is much less postweaning fallback 
in the UN/Lep cohort (Figure 2D-F). Responses also dif- 
fer greatly between adult males and females. Figure 3A 
shows the absolute weight gain on HF relative to normal 
chow diet for each group of males and females. For AD/ 
Sal, both males and females gain comparable amounts 
on a HF diet (37 g and 43 g respectively). Similarly, for 
UN/Sal, both males and females gain a larger amount on 
the HF diet (63 g and 70 g respectively). For AD/Lep, 
females gain 32 g on the HF diet while males gain 
100 g. Finally, for UN/Lep, females show an efficient res- 
cue, bringing the weight gain back down to 35 g (from 63 
in UN/Sal): by contrast males show no rescue, instead 
gaining 109 g, the highest of any treatment cohort. 
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Figure 3 Comparison of selected male and female growth data. 

A: Bar chart showing average weight in g +/- s.e.m. for all cohorts 
at day 108. Numbers in grey indicate the absolute weight gain on 
HF diet relative to chow diet for the AD/Sal, AD/Lep, UN/Sal and 
UN/Lep groups. B: Bar chart showing weight of UN/Sal/Chow (U) 
relative to AD/Sal/Chow (A) at day 30 and day 1 08. Left = females, 
right = males. 



The pathogenic programming effects of maternal un- 
dernutrition are also more severe in males than females, 
in that they are not restricted to males on the HF diet. 
Figure 3C shows the weight of UN/Sal/Chow cohorts 
relative to control AD/Sal/Chow for males and females at 
day 30 and at the end of the experimental period. Females 
and males both showed a substantial growth restriction 
(see also Figure 2C, F) that is somewhat more severe in 
males (82.9% ±3.7%) than females (88.0% ± 1.3%) at day 
30. These males then showed a continued slow weight 
gain throughout life relative to control, reaching 96.1 ± 
4.8% of control by day 108-note that this was the end of 
the male experimental series. In contrast, the females 
maintained a stable weight relative to control from day 30 
onwards (92.3% ± 1.3% on day 108, and 89.8% ± 2.4% on 
day 171-the end of the female experiment). 

Collectively, these results suggest that while the patho- 
genic "thrift- inducing" effects of early growth restriction 
are shared between males and females, the non-hepatic 
mechanisms underlying the leptin rescue in UN/Lep co- 
horts may be less effective (or even absent) in males com- 
pared to females. In the Discussion (below) we propose a 
revised and extended model for developmental program- 
ming, which reconciles the hitherto-conflicting male and 
female growth data. 

Genes associated with the prevention of metabolic 
syndrome by leptin treatment in the UN/Lep/HF cohort 

Although the AD/Lep cohorts globally show many of 
the same metabolic and transcriptional changes as the 
UN/Sal cohorts, nevertheless in females metabolic syn- 
drome is only triggered in UN/Sal/HF cohort and is pre- 
vented by postnatal leptin treatment ("leptin rescue") in 
UN/Lep/HF. It is important to determine the hepatic 
transcriptional changes associated with the leptin rescue, 
since these may be useful biomarkers for metabolic syn- 
drome and/or for the predisposed state seen in UN/Sal/ 
Chow. Genes of interest in this regard are those that show 
a significant perturbation in expression associated with 
metabolic syndrome, which is exacerbated by maternal 
undernutrition and rescued by postnatal leptin treatment. 
To identify these genes, we took the initial set of 2221 
genes showing at least some significant change in the one- 
way ANOVA and imposed three further filters: 

• Maternal diet and leptin treatment both called as 
significant, either individually or as part of an 
interaction term. 1010 genes passed this filter; the 949 
genes with a significant AB interaction plus a further 61 
genes where both factors were individually significant 
but the AB interaction term was not significant. 

• At least a 1.25 fold change in expression level between 
the affected UN/Sal/HF cohort and the control 
AD/Sal/Chow cohort. 177 genes passed this filter. 
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• Gene expression levels in the UN/Sal/HF cohort 
with metabolic syndrome must fall at an extreme 
(highest or lowest expression amongst all 8 
experimental cohorts). This yielded a final set of 98 
genes, 56 of which were upregulated in metabolic 
syndrome and 42 downregulated (Table 3 and 
Additional file 5: Table S5). 

DAVID analysis of the resulting gene list (see Additional 
file 5: Table S5) reveals significant deregulation of carbo- 
hydrate metabolic pathways (enrichment factor 2.939, 
p-value = 0.00115), steroid response genes (enrichment 
factor 1.59, p-value = 0.0256) and mitochondrially- targeted 
genes (enrichment factor 1.53, p-value = 0.0298). 

There was no significant enrichment for functional 
terms relating to immune function in this list since the 
significant downregulation of MHC class I was more pro- 
nounced in AD/Lep cohorts than in UN/ Sal cohorts and 
hence these were excluded by the final filter (the same ap- 
plies to the upregulation of complement gene C2). It is at 
present unclear why postnatal leptin treatment should 
have such a marked effect on immune parameters. 

Precursor changes in chow-fed cohorts may be associated 
with disease predisposition 

Of the four chow-fed cohorts, none directly shows meta- 
bolic syndrome, however the UN/ Sal/ Chow cohort is 
predisposed to metabolic syndrome and will suffer from 
it if fed a HF diet. We therefore looked specifically 
within the chow-fed cohorts to determine whether any 
of the 98 genes associated with metabolic syndrome 
(and rescued by leptin treatment) showed precursor ex- 
pression changes that might be involved in susceptibility 
to or initiation of metabolic syndrome. In all, 41/98 of 
these genes also showed at least a 1.25 fold expression 
change in the UN/ Sal/ Chow cohort relative to the con- 
trol AD/Sal/Chow cohort (highlighted in bold in Table 3). 
For these precursor changes, DAVID analysis did not 
show any enrichment for mitochondrially-targeted genes 
or for steroid response genes, however the annotation 
group relating to carbohydrate metabolism remained sig- 
nificant (enrichment factor 2.71, p-value = 0.002). 

Discussion 

In this study, we present a whole-genome fully factorial 
examination of the transcriptional consequences of the 
interactions between maternal diet, leptin treatment, 
postweaning diet and metabolic syndrome in a rat model 
of developmental programming. 

Our data confirm the prior observation of "leptin re- 
versal", in that of the genes regulated by postnatal leptin 
treatment, the vast majority show the opposite direction 
of leptin-induced change in pups born to undernour- 
ished mothers, when compared to those born to normally- 



nourished mothers. Strikingly, it is not simply that 
maternal undernutrition reverses the direction of the 
leptin response: rather it appears that both maternal 
undernutrition and postnatal leptin treatment inde- 
pendently switch the liver into an altered state. Para- 
doxically, the combination of both factors abolishes this 
effect and switches the liver back into a state more closely 
resembling the control (untreated) cohorts. This pattern is 
seen in both the transcriptional data and in key pheno- 
typic parameters relating to metabolic status. 

A revised model for developmental programming of 
obesity 

Based on the transcriptional findings reported above and 
on our re-examination of existing growth data for both 
male and female offspring, we propose a revised model 
of developmental programming by maternal undernutri- 
tion and postnatal leptin treatment. In this new model, 
there are two competing programmes that can be estab- 
lished quasi-independently. The first programme is a 
thrifty hepatic metabolic programme induced by early 
life growth restriction, governing the efficiency with 
which fuel is used by the body. Importantly, in the AD/ 
Lep cohorts the thrifty programming is invoked even in 
the presence of a supraphysiological leptin stimulus, sug- 
gesting that leptin itself is not the predominant signal in- 
volved in programming of hepatic metabolic efficiency. The 
second programme is a homeostatic set point programme 
governing body composition that is established in response 
to leptin stimulation during early neonatal life. The calibra- 
tion of the set point (either lean or normal body compos- 
ition) is determined by the body's composition at the time 
of leptin stimulation, while the robustness with which 
the programme is "locked in" depends on the level of the 
stimulus. Given the known effects of leptin on appetite and 
locomotor activity, it is likely that this second programme 
is mediated at least in part by centrally controlled adjust- 
ment of calorie intake and expenditure, however we cannot 
rule out effects of leptin in other non-hepatic tissues. 

Our model is based on the emerging view that during 
development leptin functions as an anti-starvation sig- 
nal, as distinct from its adult role as a fatness and sa- 
tiety signal [23]. During early development, in normally 
nourished animals, leptin levels are high (known as the 
"leptin surge") and largely independent of fat mass. As 
development proceeds, leptin levels drop to a lower 
level governed dually by overall fat mass and satiety 
[24-26]. Early life undernutrition delays or abolishes 
the leptin surge [27,28], and nutritional recovery re- 
stores it [29] -in the latter study, higher leptin levels 
were observed in "recovered" infants than in control 
infants, which we argue may represent a delayed de- 
velopmental surge rather than a simple restoration of 
leptin secretion following fat gain. 
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Table 3 Genes which show a significant interaction between maternal nutritional status and postnatal leptin treatment 
effects, an expression change of at least 1.25 fold in the UN/Sal/HF cohort (with metabolic syndrome), and where this 
cohort is an outlier from the other experimental groups 



Symbol 



Fold change 
in USC 



Fold change 
in USH 



Significant Entrez ID 
factors 



Definition 



Ifitl 

Ppp1r3b 

Ca3 
Stac3 

TUBB2A 

LEAP2 

Slc34a2 

Pgd 

Afp 

VNN3 

Avprla 

Laptm4b 

Igf2bp3 

Aldhlbl 

Amyla 

Aadac 

Pipox 
Sez6 

Olr59 

RGD1563713 
Gpd2 

Bche 
Otc 

Dnmt3b 
Clpx 

Slc7a7 

Dpys 



2.95 

1.46 

1.29 
1.27 

1.79 

1.52 

1.45 

1.24 

1.33 

1.17 

1.42 

0.94 
1.23 
1.08 

1.33 

1.19 

1.56 
1.28 

1.31 
1.31 

1.28 

1.20 
1.32 

0.95 
1.35 

1.03 

1.04 



3.52 

2.40 

2.38 
2.35 

2.33 

2.29 

1.88 

1.77 

1.76 

1.70 

1.70 

1.68 
1.67 
1.64 

1.64 

1.62 

1.60 
1.58 

1.56 
1.56 

1.56 

1.51 
1.51 

1.51 
1.50 

1.49 

1.49 



AB 56824 Rattus norvegicus interferon-induced protein with tetratri- 

copeptide repeats 1 (Ifitl ), mRNA. 

AB, C 192280 Rattus norvegicus protein phosphatase 1, regulatory 

(inhibitor) subunit 3B (Ppp1r3b), mRNA. 

AB, C 54232 Rattus norvegicus carbonic anhydrase 3 (Ca3), mRNA. 

AB, C 362895 PREDICTED: Rattus norvegicus SH3 and cysteine rich 

domain 3 (predicted) (Stac3_predicted), mRNA. 

AB, BC 498736 PREDICTED: Rattus norvegicus similar to tubulin, beta 2 

(TUBB2A), mRNA. 

AB, C 497901 PREDICTED: Rattus norvegicus similar to Liver-expressed 

antimicrobial peptide 2 precursor (LEAP-2) (LEAP2), 
mRNA. 

AB, C 84395 Rattus norvegicus solute carrier family 34 (sodium 

phosphate), member 2 (Slc34a2), mRNA. 

AB 362660 PREDICTED: Rattus norvegicus phosphogluconate 

dehydrogenase (Pgd), mRNA. 

AB, C 24177 Rattus norvegicus alpha-fetoprotein (Afp), mRNA. 

AB, C 498992 PREDICTED: Rattus norvegicus similar to Vanin-3 (predicted) 

(VNN3_predicted), mRNA. 

AB, C 25107 Rattus norvegicus arginine vasopressin receptor 1 A 

(Avprla), mRNA. 

AB 315047 Rattus norvegicus lysosomal-associated protein transmembrane 

4B (Laptm4b), mRNA. 

AB, C 312320 PREDICTED: Rattus norvegicus insulin-like growth factor 2, bind- 

ing protein 3 (Igf2bp3), mRNA. 

AB, C 298079 Rattus norvegicus aldehyde dehydrogenase 1 family, member 

B1 (Aldhlbl), nuclear gene encoding mitochondrial protein, 
mRNA. 

AB, C 24203 Rattus norvegicus amylase, alpha 1 A (salivary) (Amyla), 

mRNA. 

AB, C 57300 Rattus norvegicus arylacetamide deacetylase (esterase) (Aadac), 

mRNA. 

AB, BC 303272 Rattus norvegicus pipecolic acid oxidase (Pipox), mRNA. 

AB, C 192247 PREDICTED: Rattus norvegicus seizure related 6 homolog 

(mouse) (Sez6), mRNA. 

AB, C 170816 Rattus norvegicus olfactory receptor 59 (Olr59), mRNA. 

AB 316550 PREDICTED: Rattus norvegicus similar to Rab18 

(RGD1563713), mRNA. 

AB, C 25062 Rattus norvegicus glycerol-3-phosphate dehydrogenase 2, 

mitochondrial (Gpd2), nuclear gene encoding mitochondrial 
protein, mRNA. 

AB, C 65036 Rattus norvegicus butyrylcholinesterase (Bche), mRNA. 

AB, C 2561 1 Rattus norvegicus ornithine carbamoyltransferase (Otc), 

nuclear gene encoding mitochondrial protein, mRNA. 

AB, AC, BC 444985 Rattus norvegicus DNA methyltransferase 3B (Dnmt3b), mRNA. 

AB, BC 300786 Rattus norvegicus CIpX caseinolytic peptidase X homolog 

(E. coli) (Clpx), mRNA. 

AB, C 83509 Rattus norvegicus solute carrier family 7 (cationic amino acid 

transporter, y + system), member 7 (Slc7a7), mRNA. 

AB, C 65135 Rattus norvegicus dihydropyrimidinase (Dpys), mRNA. 



Ellis et al. BMC Genomics 2014, 15:49 
http://www.biomedcentral.eom/1 471 -21 64/1 5/49 



Page 11 of 18 



Table 3 Genes which show a significant interaction between maternal nutritional status and postnatal leptin treatment 
effects, an expression change of at least 1.25 fold in the UN/Sal/HF cohort (with metabolic syndrome), and where this 
cohort is an outlier from the other experimental groups (Continued) 



Fam82a 

Slc16a1 

Adhfel 

Inhba 
TMEM70 

RGD1 560797 

RGD1310209 

Maob 

Lrp5 

Klrgl 

Atf4 

Phfll 

Rora 

Serpina3m 

PLBD1 

ANGPTL3 

Agtrla 

Apon 
Lgals5 

Cabd 

Sdcbp 

CP 

Irf6 

Hsd17b13 
Mstl 
Btg1 
MUT 



1.34 

1.15 

1.32 

1.10 
1.13 

1.08 

1.08 

1.24 

1.09 

1.10 

1.24 

1.15 

1.25 

1.17 

1.15 

1.17 

1.02 

1.09 
1.21 

1.23 

1.12 
1.07 
1.14 

1.07 

1.17 

1.13 

1.19 



1.48 

1.41 

1.40 

1.40 
1.38 

1.38 

1.35 

1.34 

1.34 

1.34 

1.33 

1.33 

1.33 

1.32 

1.31 

1.31 

1.30 

1.30 
1.29 

1.28 

1.27 
1.27 
1.26 

1.26 

1.26 

1.25 

1.25 



AB, C 313840 Rattus norvegicus family with sequence similarity 82, 

member A (Fam82a), mRNA. 

AB 25027 Rattus norvegicus solute carrier family 16 (monocarboxylic acid 

transporters), member 1 (Slc16a1), mRNA. 

AB, C 362474 Rattus norvegicus alcohol dehydrogenase, iron 

containing, 1 (Adhfel), mRNA. 

AB, BC 29200 Rattus norvegicus inhibin beta-A (Inhba), mRNA. 

AB, C 500384 PREDICTED: Rattus norvegicus similar to RIKEN cDNA 

2210416 J16 (predicted) (TMEM70_predicted), mRNA. 

A, BC 3061 15 PREDICTED: Rattus norvegicus similar to glyceraldehyde-3- 

phosphate dehydrogenase (predicted) (RGD1560797_pre- 
dicted), mRNA. 

AB, C 36201 9 PREDICTED: Rattus norvegicus similar to KIAA1 324 protein 

(predicted) (RGD1310209_predicted), mRNA. 

AB 25750 Rattus norvegicus monoamine oxidase B (Maob), nuclear gene 

encoding mitochondrial protein, mRNA. 

AB 293649 PREDICTED: Rattus norvegicus low density lipoprotein receptor- 

related protein 5 (predicted) (Lrp5_predicted), mRNA. 

AB, C 58975 Rattus norvegicus killer cell lectin-like receptor subfamily G, 

member 1 (Klrgl), mRNA. 

AB 79255 Rattus norvegicus activating transcription factor 4 (tax- 

responsive enhancer element B67) (Atf4), mRNA. 

AB, C 361051 PREDICTED: Rattus norvegicus PHD finger protein 1 1 

(predicted) (Phfl 1_predicted), mRNA. 

AB 300807 PREDICTED: Rattus norvegicus RAR-related orphan receptor 

alpha (predicted) (Rora_predicted), mRNA. 

AB, C 299276 PREDICTED: Rattus norvegicus serine (or cysteine) proteinase 

inhibitor, clade A, member 3 M (Serpina3m), mRNA. 

AB, C 297694 Rattus norvegicus similar to RIKEN cDNA 1 100001 H23 (PLBD1), 

mRNA. 

AB, C 502970 PREDICTED: Rattus norvegicus similar to angiopoietin-related 

protein 3 (ANGPTL3), mRNA. 

AB 24180 Rattus norvegicus angiotensin II receptor, type 1 (AT1A) 

(Agtrla), mRNA. 

AB, C 304603 Rattus norvegicus apolipoprotein N (Apon), mRNA. 

AB 25475 Rattus norvegicus lectin, galactose binding, soluble 5 (Lgals5), 

mRNA. 

AB 360887 Rattus norvegicus chaperone, ABC1 activity of bcl complex 

homolog (S. pombe) (Cabd), nuclear gene encoding 
mitochondrial protein, mRNA. 

AB 83841 Rattus norvegicus syndecan binding protein (Sdcbp), mRNA. 

AB 294942 PREDICTED: Rattus norvegicus hypothetical CP (CP), mRNA. 

AB 364081 PREDICTED: Rattus norvegicus interferon regulatory factor 6 

(predicted) (Irf6_predicted), mRNA. 

AB, C 305150 PREDICTED: Rattus norvegicus hydroxysteroid (17-beta) 

dehydrogenase 13 (Hsd17b13), mRNA. 

AB 24566 Rattus norvegicus Macrophage stimulating 1 (hepatocyte 

growth factor-like) (Mstl), mRNA. 

AB, C 29618 Rattus norvegicus B-cell translocation gene 1, anti-proliferative 

(Btg1), mRNA. 

AB 497857 PREDICTED: Rattus norvegicus similar to MYLE protein 

(Dexamethasone-induced protein) (predicted) (Mispredicted), 
mRNA. 
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Table 3 Genes which show a significant interaction between maternal nutritional status and postnatal leptin treatment 
effects, an expression change of at least 1.25 fold in the UN/Sal/HF cohort (with metabolic syndrome), and where this 
cohort is an outlier from the other experimental groups (Continued) 



Pxmp4 
Igfbp2 
Sds 

Tnfsf13 

Hnrpab 

Asl 
Cd63 

RGD1 560687 
Trpm6 

Hspbl 

Ftl1 

Mgst2 

Znf593 

Ppcs 

Mthfdl 

Pmml 

Aqpll 
Csrp2 

Pla2g12a 

Tifa 

Hexb 
Rps14 
Ier3 

Paics 

Txnl5 
Plala 

Cyp8b1 



-1.38 

-2.11 

-1.54 
-2.22 

-1.19 

-1.33 
-1.36 
-1.25 

-1.27 

-1.34 

-1.11 
-1.43 

-1.19 

-1.43 

-1.38 

-1.26 

-1.35 
-1.34 

-1.27 

-1.24 

-1.17 
-1.15 
-1.38 

-1.40 

-1.32 
-1.35 

-1.08 



-3.51 

-3.09 

-2.49 
-2.42 

-2.09 

-2.07 
-1.94 
-1.85 

-1.76 

-1.68 

-1.60 
-1.54 

-1.53 

-1.52 

-1.47 

-1.47 

-1.45 
-1.45 

-1.44 

-1.42 

-1.41 
-1.41 
-1.40 

-1.40 

-1.38 
-1.38 

-1.37 



AB 282634 Rattus norvegicus peroxisomal membrane protein 4 

(Pxmp4), mRNA. 

AB, C 25662 Rattus norvegicus insulin-like growth factor binding 

protein 2 (Igfbp2), mRNA. 

AB 25044 Rattus norvegicus serine dehydratase (Sds), mRNA. 

AB 287437 Rattus norvegicus tumor necrosis factor (ligand) 

superfamily, member 13 (Tnfsf13), mRNA. 

AB, C 83498 Rattus norvegicus heterogeneous nuclear ribonucleoprotein 

A/B (Hnrpab), mRNA. 

AB, C 59085 Rattus norvegicus argininosuccinate lyase (Asl), mRNA. 

AB, C 29186 Rattus norvegicus CD63 antigen (Cd63), mRNA. 

AB, C 500804 PREDICTED: Rattus norvegicus similar to Ferritin light 

chain (Ferritin L subunit) (predicted) 
(RGD1560687_predicted), mRNA. 

AB, C 293874 PREDICTED: Rattus norvegicus transient receptor potential 

cation channel, subfamily M, member 6 (predicted) 
(Trpm6_predicted), mRNA. 

AB, C 24471 Rattus norvegicus heat shock 27 kDa protein 1 (Hspbl), 

mRNA. 

AB, C 29292 Rattus norvegicus ferritin light chain 1 (Ftll), mRNA. 

AB, BC 295037 PREDICTED: Rattus norvegicus microsomal glutathione 

S-transferase 2 (predicted) (Mgst2_predicted), mRNA. 

AB, C 298546 PREDICTED: Rattus norvegicus zinc finger protein 593 

(predicted) (Znf593_predicted), mRNA. 

AB, C 298490 Rattus norvegicus phosphopantothenoylcysteine 

synthetase (Ppcs), mRNA. 

A, BC 64300 Rattus norvegicus methylenetetrahydrofolate 

dehydrogenase (NADP + dependent) 1, 
methenyltetrahydrofolate cyclohydrolase, 
formyltetrahydrofolate synthetase (Mthfdl), mRNA. 

A, B, C 300089 Rattus norvegicus phosphomannomutase 1 (Pmml), 

mRNA. 

AB 286758 Rattus norvegicus aquaporin 11 (Aqpll), mRNA. 

A, B, C 29317 Rattus norvegicus cysteine and glycine-rich protein 2 

(Csrp2), mRNA. 

AB, BC 362039 PREDICTED: Rattus norvegicus phospholipase A2, group 

XIIA (predicted) (Pla2g12a_predicted), mRNA. 

AB, C 310877 Rattus norvegicus TRAF-interacting protein with 

forkhead-associated domain (Tifa), mRNA. 

AB, C 294673 Rattus norvegicus hexosaminidase B (Hexb), mRNA. 

AB, C 29284 Rattus norvegicus ribosomal protein S14 (Rps14), mRNA. 

AB 294235 Rattus norvegicus immediate early response 3 (Ier3), 

mRNA. 

AB 140946 Rattus norvegicus phosphoribosylaminoimidazole 

carboxylase, phosphoribosylaminoimidazole 
succinocarboxamide synthetase (Paics), mRNA. 

AB 287474 PREDICTED: Rattus norvegicus thioredoxin-like 5 

(predicted) (Txnl5_predicted), mRNA. 

AB 85311 Rattus norvegicus phospholipase Al member A (Plala), 

mRNA. 

AB, C 81924 Rattus norvegicus cytochrome P450, family 8, subfamily b, 

polypeptide 1 (Cyp8b1), mRNA. 
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Table 3 Genes which show a significant interaction between maternal nutritional status and postnatal leptin treatment 
effects, an expression change of at least 1.25 fold in the UN/Sal/HF cohort (with metabolic syndrome), and where this 
cohort is an outlier from the other experimental groups (Continued) 



Creb3 


-1.28 


-1.36 


AB, C 


298400 


Rattus norvegicus cAMP responsive element binding 
protein 3 (Creb3), mRNA. 


Cd320 


-1.06 


-1.36 


AB, C 


362851 


PREDICTED: Rattus norvegicus CD320 antigen (Cd320), mRNA. 


Nola2 


-1.10 


-1.34 


AB, C 


287273 


PREDICTED: Rattus norvegicus nucleolar protein family A, 
member 2 (predicted) (Nola2_predicted), mRNA. 


Sultlcl 


-0.87 


-1.32 


AB 


65185 


Rattus norvegicus sulfotransferase family, cytosolic, 1C, 
member 1 (Sultlcl), mRNA. 


Slc11a2 


-1.16 


-1.32 


AB, BC 


25715 


Rattus norvegicus solute carrier family 1 1 (proton-coupled 
divalent metal ion transporters), member 2 (Slcl 1a2), mRNA. 


Btbd9 


-1.19 


-1.31 


A, BC 


294318 


Rattus norvegicus BTB (POZ) domain containing 9 (Btbd9), 
mRNA. 


Htatip2 


-1.20 


-1.30 


AB, C 


292935 


PREDICTED: Rattus norvegicus HIV-1 tat interactive protein 2, 
homolog (human) (predicted) (Htatip2_predicted), mRNA. 


CtsM 


-1.09 


-1.29 


AB, C 


25697 


Rattus norvegicus cathepsin L1 (CtsM), mRNA. 


Psmb3 


-1.14 


-1.29 


AB, C 


29676 


Rattus norvegicus proteasome (prosome, macropain) subunit, 
beta type 3 (Psmb3), mRNA. 


Sat2 


-1.16 


-1.29 


AB 


360547 


PREDICTED: Rattus norvegicus spermidine/spermine N1 -acetyl 
transferase 2 (predicted) (Sat2_predicted), mRNA. 


Gtf2b 


-1.15 


-1.29 


AB, C 


81673 


Rattus norvegicus general transcription factor IIB (Gtf2b), 
mRNA. 


Eppb9 


-1.10 


-1.26 


AB, BC 


287383 


PREDICTED: Rattus norvegicus endothelial precursor protein B9 
(predicted) (Eppb9_predicted), mRNA. 


Sdf2 


-1.05 


-1.26 


AB 


287470 


PREDICTED: Rattus norvegicus stromal cell derived factor 2 
(predicted) (Sdf2_predicted), mRNA. 


Wdr45 


-1.17 


-1.26 


AB, C 


302559 


Rattus norvegicus WD repeat domain 45 (Wdr45), mRNA. 


Psma5 


-1.20 


-1.26 


AB, C 


29672 


Rattus norvegicus proteasome (prosome, macropain) subunit, 
alpha type 5 (Psma5), mRNA. 



Bold indicates genes where the expression change in the UN/Sal/Chow cohort (predisposed to metabolic syndrome) is also at least 1.25 fold, indicating that 
precursor changes are occurring in this group also. Additional file 5: Table S5 is an expanded version of this table containing further annotation relating to gene 
ontology, gene function and known human diseases related to these genes. 



We hypothesise that the function of neonatal leptin 
signalling is to train the homeostatic set point mecha- 
nisms to recognise and subsequently maintain an appro- 
priate body fat composition, and that this pattern (surge 
followed by drop) represents the calibration mechanism 
for the feedback loop. The peak of the leptin surge trig- 
gers developments underpinning set point programming, 
while the leptin levels immediately subsequent to the 
surge set the baseline which will subsequently be 
defended by the newly- calibrated set point program. If 
the individual is normal weight at the time of the leptin 
surge, a normal calibration is established and will subse- 
quently be maintained. Under starvation conditions the 
calibration process would normally be delayed by mech- 
anisms including lower overall leptin secretion and also 
competition by excess soluble leptin receptor [29] -how- 
ever if exogenous leptin is administered, the combination 
of a leptin surge with an underweight body composition 
means that an unusually lean calibration is established 
and will subsequently be maintained. Conversely, if the 



individual is overweight at the time of the leptin surge, the 
calibration is set inappropriately high. 

Interpreting the differences between male and female 
responses to early programming 

Our existing growth data for male rats suggests that the 
homeostatic set point mechanism is less efficient in 
males than in females. Interestingly, in human babies 
there is a sex difference in favour of higher leptin con- 
centrations in female newborns [26], which is consistent 
with our above hypothesis about the function of the lep- 
tin surge in set point calibration. Incorporating this final 
factor leads to the unified model described in Table 4, 
which summarises the key elements of our new model 
and how it relates to the observed phenotypes in both 
sexes. Briefly, the first, hepatic element of the programme 
(i.e. thrifty use of fuel) is induced by early life growth 
restriction irrespective of the reason for that growth re- 
striction, and as such all except the AD/Sal cohorts are 
programmed to be thrifty. The second, leptin-mediated 
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Table 4 Unified model for how the competing hepatic thriftiness and non-hepatic set point programmes interact to 
produce the observed phenotypes in males and females 



Cohort Metabolism 1 
$ AD/Sal Normal 
AD/Lep Thrifty 



UN/Sal Thrifty 

UN/Lep Thrifty 

S AD/Sal Normal 

AD/Lep Thrifty 

UN/Sal Thrifty 

UN/Lep Thrifty 



Set point calibration 

Normal 

Normal 



Lean 

Lean 

Normal 
Normal 

Lean 

Lean 



Set point enforcement 3 Phenotypic consequences 

Strong (female default) Baseline for chow and high fat diet regimes in females. 



Strong (female default) 



Impaired (due to 
insufficient leptin 
availability) 

Strong (restored to 
female default level) 



Weak (male default) 
Weak (male default) 



Severely impaired 

(weak male default is 
lowered further) 

Weak (restored to 
male default level) 



Normal set point calibration is enforced and overrides the thrifty 
metabolism. Weight remains comparable to AD/Sal females on 
both low and high fat diets despite thrifty liver biochemistry. 

Set point enforcement is impaired and consequently unable to 
compensate for a thrifty metabolism when fed a HF diet. UN/Sal/HF 
females become morbidly obese and suffer metabolic syndrome. 

A well-enforced "lean" set point in rescued females compensates for the 
thrifty metabolism, enforcing a lean body composition in UN/Lep/Chow 
and keeping weight gain down to normal levels in UN/Lep/HF. 

Baseline for chow and high fat diet regimes in males. 

As with UN/Sal animals (male and female), set point enforcement 
is weak and cannot compensate for the thrifty metabolism. 
The AD/Lep HF phenotype is more severe than UN/Sal/HF 
because the set point is higher. 

Set point enforcement is severely impaired and thus UN/Sal males 
cannot maintain a steady body composition even on a normal 
chow diet. However, the set point is still lean and so the 
phenotype on HF diet is less severe than AD/Lep males. 

Combined programme is similar to UN/Sal females. Rescued 
males maintain a steady body composition on a normal chow 
diet but show elevated weight gain on HF diet. 



The thrifty metabolic programme is induced by periods of restricted growth, whether due to maternal undernutrition or postnatal leptin treatment, and is 
consequently present in all except AD/Sal cohorts. 

2 The level of the set point is established by the body composition at the time of the neonatal leptin surge, and is consequently normal for AD cohorts and lean 
for UN cohorts. 

3 The strength with which the set point calibration is enforced depends on the levels of leptin at the time of calibration. Set point enforcement is weaker in males 
than in females for unknown reasons. 



element of the programme (i.e. the homeostatic set point) 
varies considerably between cohorts. The AD/Sal and AD/ 
Lep cohorts have a normal body composition at the time 
of leptin stimulation and thus are calibrated to a normal 
set point. The UN/Lep cohorts have a very lean body 
composition at the time of leptin stimulation and thus are 
calibrated to a lean set point. The UN/Sal cohorts also 
have a lean body composition postnatally and are there- 
fore calibrated to a lean set point, but since there is no ex- 
ogenous leptin given and endogenous levels are low, the 
calibration is only weakly enforced. Consequently, in the 
UN/Lep cohort the set point mechanisms are strong 
enough to over-ride the metabolic programme and 
enforce a lean body composition, masking the thrifty 
hepatic metabolic programming, while in the UN/ Sal co- 
hort they cannot do so. In males the hepatic element of 
the programme is maintained, while the set point enforce- 
ment is weaker for unknown reasons. 

The new model thus explains four key puzzles con- 
cerning various features of the programmed phenotypes 
in males and females of this rat model, namely; 

(1) In females AD/Lep and UN/Sal both induce 

transcriptional changes related to fatty liver disease, 
but the combination of both in UN/Lep does not 
(Figure 1 and Results passim). 



(2) AD/Lep treatment is obesogenic in males but not in 
females (Figure 3A). 

(3) Neonatal leptin treatment protects against metabolic 
syndrome in UN females on a high fat postweaning 
diet, but is much less effective in males (Figure 3A). 

(4) UN/Sal/Chow females maintain a constant weight 
relative to AD/Sal/Chow controls from weaning 
onwards, whereas UN/Sal/Chow males shows a 
lifelong slow increase relative to AD/Sal/Chow 
controls (Figure 3B). 

Significance of the "calibration" model for postnatal 
leptin signalling 

Potentially this view of postnatal leptin signalling as a 
training/calibration mechanism is also applicable in wider 
contexts. For example, maternal overnutrition (as well as 
undernutrition) can also programme metabolic syndrome 
and obesity in offspring [21]. Under our model, this is ex- 
plained by the offspring being fatter at the time of the 
postnatal leptin training signal. The set point is conse- 
quently established at a higher level, leading these individ- 
uals to maintain and defend an inappropriately high body 
fat composition. 

Similarly, this model could explain the differences in 
outcome (seen in both human and animal models) be- 
tween neonates exhibiting rapid catch-up growth after 



Ellis et al. BMC Genomics 2014, 15:49 
http://www.biomedcentral.eom/1 471 -21 64/1 5/49 



Page 15 of 18 



starvation and those with a slower return to normal 
weight [30,31]. Given that rapid weight recovery leads 
disproportionately to fat gain rather than lean tissue 
gain, it will skew the overall body composition to a higher 
fat percentage, which will subsequently be "locked in" 
once the leptin signalling reaches a high enough level to 
trigger the establishment of a body composition set point. 
Consistent with this, in a rat model, rapid catch up growth 
leads to leptin resistance in adulthood-the reduced sensi- 
tivity will consequently lead these animals to self-regulate 
to a higher plasma leptin level and a higher body fat con- 
tent. Conversely, slower catch up growth does not have 
the same effect on the leptin feedback loop [32,33]. 

Of key interest for future research will be whether this 
post-natal calibration of the body's set point is a one-off 
event during a specific developmental window, or whether 
it can be re-triggered by a leptin surge (endogenous or ex- 
ogenous) at later ages. One can envisage a "ratchet" mech- 
anism for morbid obesity whereby weight gain beyond a 
given threshold triggers a re- calibration and establishment 
of a "new normal", making subsequent weight loss harder. 
Conversely, if "surge" levels of leptin were administered 
following a period of weight loss, would it then be able to 
re-establish a healthier set point for the patient? 

A thrifty liver is an immunosuppressed, inflamed liver 

In functional terms, there are two core axes affected by 
the phenotypic switch in the liver, namely metabolism- 
and immune-related functions. The metabolism-related 
changes include upregulation of multiple genes concerned 
with mitochondrial function, glycolysis, gluconeogenesis 
and carbohydrate metabolism. This is associated with in- 
creased total body fat percentage, increased circulating 
fasting leptin levels, and increased insulin/C-peptide levels. 

The second functional cluster of expression changes 
observed relates to immune function, with a concerted 
downregulation of MHC genes (predominantly class I 
genes, but also some class lb and class II genes) and the 
Tap2 peptide transporter necessary for antigen presenta- 
tion by class I molecules. Concurrent with this downreg- 
ulation of adaptive immunity, there is a trend towards 
upregulation of the innate immune system, including both 
generic inflammatory mediators such as complement and 
also liver-specific markers such as alpha-fetoprotein. These 
changes are consistent with a model of non-alcoholic 
fatty liver disease as a chronic low-grade inflammatory 
process [34]. 

Concurrent with the pro-inflammatory changes ob- 
served, the effects on the adaptive immune system are 
interesting in relation to the immunological effects of 
perinatal programming, in particular the increased mor- 
tality from infectious disease in in utero-depnved chil- 
dren [13,14]. Mechanistically, there is a plausible link 
between the metabolic phenotype and the immune system 



changes in that excess saturated fatty acids are known to 
exert an immunosuppressive effect via reduced antigen 
presentation by MHC class I molecules [35-37]. Poten- 
tially therefore this rat model may in future serve as a 
model for immune system development in malnourished 
children, as well as for metabolic disease. 

It is intriguing to note that these transcriptional 
changes-both those reflecting increased inflammation 
and those reflecting decreased antigen presentation-are 
not only seen in the UN/Sal/HF cohort, but are consist- 
ent across all AD/Lep and UN/Sal cohorts irrespective 
of postweaning diet. This suggests that the immune 
changes are directly associated with the thrifty pheno- 
type itself rather than with metabolic syndrome. An im- 
portant avenue for future research will therefore be to 
examine immune parameters at early timepoints in pro- 
grammed thrifty and non-thrifty animals to confirm 
whether the immune changes do indeed precede the de- 
velopment of clinical disease, and what impact this has 
for our understanding of the pathogenesis of obesity. 

Do variations in hepatic thriftiness affect predisposition 
to metabolic syndrome in a nutritionally rich 
environment? 

In the female experiment analysed here, despite the fact 
that both AD/Lep and UN/Sal provoke the same thrifti- 
ness phenotypic switch within the liver, it is only the 
UN/Sal rats that subsequently develop the full hallmarks 
of metabolic syndrome when fed a high fat diet. As we 
argue above, we believe this is most likely due to non- 
hepatic effects of leptin in establishing and maintaining 
a set point for body composition, however it is also 
plausible that there may be varying degrees of hepatic 
thriftiness depending on the nature and intensity of the 
stimulus provoking the programmed changes. 

Therefore, we looked in detail at genes most strongly 
associated with pathological status (i.e. genes where UN/ 
Sal/HF is an outlier from all other cohorts) and which 
were also affected by the maternal diet/leptin inter- 
action. These are the genes most strongly associated 
with the induction of metabolic syndrome by maternal 
undernutrition and its subsequent rescue by leptin treat- 
ment. These too showed disturbances of carbohydrate 
metabolism and an increase in gluconeogenesis. Strik- 
ingly, around a quarter of this final set of genes showed 
precursor changes in the UN/ Sal/ Chow cohort, which 
potentially represent pre-pathogenic events that predis- 
pose to metabolic syndrome even among rats on a nor- 
mal diet. Further investigation of this cohort could 
potentially indicate biomarkers for those at risk of meta- 
bolic syndrome, allowing targeted intervention and/or 
advice to be given before the onset of pathological 
change. Conversely, investigation of the non-hepatic ef- 
fects of leptin in the AD/Lep/HF cohort (at least in 
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females) may identify the factors that protect against 
metabolic syndrome even in the presence of a thrifty 
hepatic phenotype. 

Conclusions 

In this study, we have established that both prenatal un- 
dernutrition and postnatal leptin treatment lead to a 
thrifty liver programme, and explain this as a nonspecific 
response to growth restriction from any cause. We fur- 
ther show that this thrifty programming is associated 
with increased expression of inflammatory markers and 
also with a downregulation of antigen presentation genes 
that may thus lead to immunosuppression. We propose 
a new model for developmental programming of obesity 
involving competing thriftiness and set point programmes, 
which resolves several conflicts between the data from 
male and female rats, and may apply more generally to 
other models of obesogenic programming. 

Methods 

Liver samples 

The experimental design utilised in this study has been 
described previously [17-19]. Briefly, virgin Wistar rats 
were time-mated and assigned to two nutritional groups: 
ad-libitum (AD) or undernourished (UN, 30% of ad- 
libitum). After birth, female AD and UN pups were ran- 
domized to receive either saline (Sal) or recombinant rat 
leptin (Lep, 2.5 (ig/g-d) on neonatal days 3-13. After 
weaning (day 22), saline-or leptin-treated AD and UN off- 
spring were weight matched and placed on either standard 
rat chow or a high fat (HF) chow (Research Diets No. 
12451; 45% energy as fat) for the remainder of the study. 
At postnatal day 170 rats were fasted overnight and killed 
by halothane anaesthesia followed by decapitation. This 
results in a balanced factorial 2x2x2 design with 8 treat- 
ment groups, n = 8 animals per group. Liver samples were 
immediately removed (collected from the same lobe) and 
immediately snap frozen in liquid nitrogen and stored 
at-80°C for later analysis. Time of culling was rotated 
within and between groups to avoid circadian con- 
founders, and was restricted to a narrow window be- 
tween 9-11 am, at the start of the light phase of the 
standard 12/12 dark/light cycle. 

RNA extraction 

RNA from whole liver tissue was prepared using Trizol 
reagent (Sigma) according to the manufacturers proto- 
col. Briefly, for each individual, 25 mg of frozen liver 
tissue was homogenised in 1 ml Trizol, followed by 
chloroform extraction. The aqueous supernatant con- 
taining the RNA was retained and the RNA precipi- 
tated with isopropanol. RNA pellets were washed with 
50% ethanol, air dried and resuspended in RNAse-free 
water (Milli-Q). RNA concentration and integrity was 



assayed with a Nanodrop 1000 spectrophotometer and an 
Agilent 2100 bioanalyser. 

Microarray hybridisation and scanning 

Array profiling was performed using the RatRef-12 
oligonucleotide platform (Illumina). This comprises 21,910 
probes covering the vast majority of the rat transcrip- 
tome. Probe labelling, hybridisation, washing and scan- 
ning were performed according to the manufacturer s 
protocols using the Illumina Total Prep kit (Applied 
Biosystems). Briefly, first strand cDNA was synthesised 
in a total volume of 20 \x\ with the supplied reagents. 
The complete first strand product was used for second 
strand synthesis, followed by column purification. The 
purified product was then used for in vitro transcrip- 
tion using T7 polymerase. Biotin-16-dUTP was incorpo- 
rated during this step, resulting in a biotinylated cRNA 
(complementary RNA) probe suitable for hybridisation. 
Probe integrity was verified using the Nanodrop 100 
and Agilent 2100, as for the initial RNA samples. Labelled 
cRNA (1.5 (ig) was hybridised to the array overnight at 
55°C in a total volume of 30 \A of the manufacturers hy- 
bridisation buffer, followed by post-hybridisation stringency 
washing and scanning (BeadArray Reader, Illumina) using 
the manufacturer s standard protocols. 

Data extraction and QC 

Array data was extracted and presence/absence calls 
performed using BeadStudio (Illumina), and subsequently 
normalised in Lumi using a variance-stabilising transform- 
ation. For each set of eight measurements for each gene/ 
experimental condition combination, the single measure- 
ment furthest from the mean was excluded. This has the 
effect of robustly eliminating noisy outlier measurements 
while maintaining a fully balanced experimental design. 
After outlier exclusion, an initial filter was applied to re- 
move non-expressed genes (defined as those called as ab- 
sent in > =75% of the data set). This reduced the initial set 
of 21,910 probes to 8,497 liver-expressed genes for further 
analysis. 

ANOVA analysis, gene classification and filtering 

One-way ANOVA analysis and gene list filtering was per- 
formed using Kensington Discovery Edition (Inforsense). 
Benjamini/Hochberg correction was used to control for 
multiple testing, with a 5% false discovery rate (FDR). This 
yielded a list of 2221 probes showing significant regulation 
in at least one experimental cohort. Following Zhou et al. 
2011 [38], these genes were then classified via a three-way 
ANOVA looking at all three experimental factors (A = ma- 
ternal diet, B = leptin treatment, C = postweaning diet) 
and all four potential interaction terms (AB/BC/AC/ABC 
interactions). Three-way ANOVA was performed using 
GNU Octave. 
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Additional file 2: Table S2. All genes found to show a significant 
expression change, classified by which factor(s)/and interaction(s) were 
significant. Individual tabs are given for each factor or combination of 
factors. 

Additional file 3: Table S3. Front tab: Cluster analysis of the 949 genes 
with a significant interaction between maternal diet and leptin. Further 
tabs show the annotation groups detected by DAVID as significantly 
enriched in each cluster. 

Additional file 4: Table S4. Weight data for both male and female 
experimental series at all time points from birth to sacrifice. 

Additional file 5: Table S5. Front tab: An expanded version of Table 3 
(disease-associated genes regulated by early life programming: bold 
indicates "precursor genes" with altered expression in Un/Sal/Chow), with 
additional columns relating to known disease associations documented 
in OMIM and other databases. Remaining tabs show the annotation 
groups detected by DAVID as significantly enriched among disease- 
associated genes and precursor genes. 



First, genes with a significant three-way interaction 
after FDR correction were assigned to group ABC. Next, 
genes with no significant three-way interaction, but where 
one or more of the two-way interaction terms was signifi- 
cant after FDR correction were assigned to groups AB/ 
AC/BC as appropriate. Finally, genes with significant 
effects for each individual experimental factor and no 
confounding interaction term were assigned to groups 
A/B/C as appropriate. The categories are not necessarily 
exclusive: for example a gene showing additive non- 
interacting effects of both maternal diet and leptin treat- 
ment is included in both categories A and B. However, a 
gene with a significant interaction term, i.e. where both 
factors are significant but the effects do not combine 
linearly, would instead fall into category AB. This iterative 
method of classification ensures that a gene is only in- 
cluded in a given category if there are no confounding 
higher-order interaction terms. 

DAVID analysis 

DAVID analysis was performed using the web inter- 
face available at http://david.abcc.ncifcrf.gov/summary.jsp 
using default parameters for clustering stringency and the 
Illumina rat chipset (RatRef-12_Vl_0_R4_11222119_A) as 
the reference list from which to calculate functional clus- 
ter enrichment. 

Phenotypic measurements 

Plasma leptin was measured by in-house radioimmuno- 
assay as described previously (Vickers et al, Endocrinology 
2001). Fasting plasma insulin, C-peptide and total ghrelin 
were measured using commercially available kits (rat insu- 
lin ELISA, Cat# 10-1124-10; Mercodia, Uppsala, Sweden; 
C-peptide RIA, cat# RCP-21 K; Linco Research Inc., St. 
Charles, MO; Ghrelin RIA, cat# GHRT-89HK, Linco 
Research). Fasting plasma glucose was measured using a 
YSI glucose analyzer (Model 2300; Yellow Springs Instru- 
ment Co., OH, USA). Plasma FFAs, glycerol and triglycer- 
ides were measured in-house using standard colorimetric 
assays. 

Availability of supporting data 

The expression array data set supporting the results of this 
article is available via LabOnline, doi:10.6070/H4X63JW 
https://mynotebook.labarchives.com/share/Leptin/MjMuN 
HwyOTA2NC8xOC03LlRyZWOb2RlLzMzMzM5NTUw 
Njh8NTkuNA==. Further supplementary data files are 
provided online as follows: 

Additional files 



Additional file 1: Table SI. Phenotypic parameters for the 
experimental series and 3-factor ANOVA analysis of significance. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

PDG, MHV and NAA conceived the study. PDG and MHV provided the tissue 
samples used for array analysis and carried out the measurements of 
phenotypic parameters. TJM performed sample extraction and microarray 
hybridisation. PJIE, TJM and BMS carried out the microarray analysis and 
functional clustering. SG and CAS advised on the functional interpretation of 
array results. PJIE drafted the manuscript. PJIE, BMS, NAA, MHV and TJM 
reviewed/edited the manuscript and all authors read and approved the final 
manuscript. 

Acknowledgements 

Cambridge (PJIE, BMS, TJM, CAS, NAA): This work was supported by the 
BBSRC and the Leverhulme Trust. Auckland (MHV, SG, PDG): This work was 
supported via the Health Research Council of New Zealand and Gravida: 
National Centre for Growth and Development. We thank Julien Bauer for 
assistance with microarray data processing and normalisation. 

Author details 

University of Cambridge Department of Pathology, Tennis Court Road, 
Cambridge CB2 1QP, UK. 2 UCL Cancer Institute, University College London, 
London WC1E 6BT, UK. 3 Liggins Institute and Gravida: National Centre for 
Growth and Development, University of Auckland, Auckland, New Zealand. 

Received: 25 October 2013 Accepted: 6 January 2014 
Published: 21 January 2014 

References 

1. Gluckman P, Hanson M: The developmental origins of health and disease 
and overview. In Developmental Origins of Health and Disease. Edited by 
Gluckman P, Hanson M. Cambridge: Cambridge University Press; 2006:1-5. 

2. Gluckman PD, Hanson MA, Bateson P, Beedle AS, Law CM, Bhutta ZA, 
Anokhin KV, Bougneres P, Chandak GR, Dasgupta P, Smith GD, Ellison PT, 
Forrester TE, Gilbert SF, Jablonka E, Kaplan H, Prentice AM, Simpson SJ, 
Uauy R, West-Eberhard MJ: Towards a new developmental synthesis: 
adaptive developmental plasticity and human disease. Lancet 2009, 
373:1654-1657. 

3. Godfrey KM, Gluckman PD, Hanson MA: Developmental origins of 
metabolic disease: life course and intergenerational perspectives. 
Trends Endocrinol Metab 2010, 21:199-205. 

4. Fernandez-Twinn DS, Ozanne SE: Early life nutrition and metabolic 
programming. Ann N Y Acad Sci 2010, 1212:78-96. 

5. Li M, Sloboda DM, Vickers MH: Maternal obesity and developmental 
programming of metabolic disorders in offspring: evidence from animal 
models. Exp Diabetes Res 201 1, 201 1:592408. 

6. Patti ME: Intergenerational programming of metabolic disease: evidence 
from human populations and experimental animal models. Cell Mol Life 
Sci 2013, 70:1597-1608. 



Ellis et al. BMC Genomics 201 4, 1 5:49 Page 1 8 of 1 8 

http://www.biomedcentral.eom/1 471 -21 64/1 5/49 



10. 



12. 



13. 



14. 



16. 



17. 



19. 



20. 



21. 



22. 



23. 



24. 



25. 



26. 



27. 



29. 



30. 



Ozanne SE: Metabolic programming in animals. Br Med Bull 2001, 
60:143-152. 

Ozanne SE, Lewis R, Jennings BJ, Hales CN: Early programming of weight 
gain in mice prevents the induction of obesity by a highly palatable 
diet. Clin Sci (Lond) 2004, 106:141-145. 

Gluckman PD, Hanson MA, Beedle AS, Spencer HG: Predictive adaptive 
responses in perspective. Trends Endocrinol Metob 2008, 19:109-1 10. 
Vaag AA, Grunnet LG, Arora GP, Brans C: The thrifty phenotype hypothesis 
revisited. Diabetologia 2012, 55:2085-2088. 

Hales CN, Barker DJ: The thrifty phenotype hypothesis. Br Med Bull 2001, 
60:5-20. 

Brewster DR, Greenwood BM: Seasonal variation of paediatric diseases in 
The Gambia, West Africa. Ann Trop Poediotr 1993, 13:133-146. 
Moore SE, Cole TJ, Poskitt EM, Sonko BJ, Whitehead RG, McGregor IA, 
Prentice AM: Season of birth predicts mortality in rural Gambia. Nature 1997, 
388:434. 

Moore SE, Cole TJ, Collinson AC, Poskitt EM, McGregor IA, Prentice AM: 
Prenatal or early postnatal events predict infectious deaths in young 
adulthood in rural Africa. Int J Epidemiol 1999, 28:1088-1095. 
Moore SE, Fulford AJ, Streatfield PK, Persson LA, Prentice AM: Comparative 
analysis of patterns of survival by season of birth in rural Bangladeshi 
and Gambian populations. Int J Epidemiol 2004, 33:137-143. 
Simondon KB, Elguero E, Marra A, Diallo A, Aaby P, Simondon F: Season of 
birth is not associated with risk of early adult death in rural 
Senegal. Int J Epidemiol 2004, 33:130-136. 

Vickers MH, Gluckman PD, Coveny AH, Hofman PL, Cutfield WS, Gertler A, 

Breier BH, Harris M: Neonatal leptin treatment reverses developmental 

programming. Endocrinology 2005, 146:421 1-4216. 

Vickers MH, Gluckman PD, Coveny AH, Hofman PL, Cutfield WS, Gertler A, 

Breier BH, Harris M: The effect of neonatal leptin treatment on postnatal 

weight gain in male rats is dependent on maternal nutritional status 

during pregnancy. Endocrinology 2008, 149:1906-1913. 

Gluckman PD, Lillycrop KA, Vickers MH, Pleasants AB, Phillips ES, Beedle AS, 

Burdge GC, Hanson MA: Metabolic plasticity during mammalian 

development is directionally dependent on early nutritional status. 

Proc Natl Acad Sci USA 2007, 1 04:1 2796-1 2800. 

Huang DW, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: 

paths toward the comprehensive functional analysis of large gene lists. 

Nucleic Acids Res 2009, 37:1-13. 

Howie GJ, Sloboda DM, Kamal T, Vickers MH: Maternal nutritional history 
predicts obesity in adult offspring independent of postnatal diet. 

J Physiol 2009, 587:905-915. 

Pico C, Jilkova ZM, Kus V, Palou A, Kopecky J: Perinatal programming of 
body weight control by leptin: putative roles of AMP kinase and muscle 
thermogenesis. Am J Clin Nutr 201 1, 94:1 830S-1 837S. 
Prentice AM, Moore SE, Collinson AC, O'Connell MA: Leptin and 
undernutrition. Nutr Rev 2002, 60:S56-S67. discussion S68-84, 85-87. 
Frederich RC, Lollmann B, Hamann A, Napolitano-Rosen A, Kahn BB, Lowell 
BB, Flier JS: Expression of ob mRNA and its encoded protein in rodents: 
impact of nutrition and obesity. J Clin Invest 1995, 96:1658-1663. 
Ahima RS, Prabakaran D, Flier JS: Postnatal leptin surge and regulation of 
circadian rhythm of leptin by feeding: implications for energy homeostasis 
and neuroendocrine function. J Clin Invest 1998, 101:1020-1027. 
Collinson A, Moore S, O'Connell M, Charalambos C, Prentice A: Developmental 
changes in leptin as a measure of energy status in human infants in 
a natural ecologic setting. Am J Clin Nutr 2005, 81 :488-494. 
Bautist CJ, Boeck L, Larrea F, Nathanielsz PW, Zambrano E: Effects of a 
maternal low protein isocaloric diet on milk leptin and progeny serum 
leptin concentration and appetitive behavior in the first 21 days of 
neonatal life in the rat. Pediatr Res 2008, 63:358-363. 
Delahaye F, Breton C, Risold PY, Enache M, Dutriez-Casteloot I, Laborie C, 
Lesage J, Vieau D: Maternal perinatal undernutrition drastically reduces 
postnatal leptin surge and affects the development of arcuate 
nucleus proopiomelanocortin neurons in neonatal male rat pups. 
Endocrinology 2008, 149:470-475. 

Stein K, Vasquez-Garibay E, Kratzsch J, Romero-Velarde E, Jahreis G: Influence 
of nutritional recovery on the leptin axis in severely malnourished children. 

J Clin Endocrinol Metab 2006, 91 :1 021 -1 026. 

Hales CN, Ozanne SE: The dangerous road of catch-up growth. J Physiol 
2003, 547:5-10. 



31 . Ong KK, Ahmed ML, Emmett PM, Preece MA, Dunger DB: Association 
between postnatal catch-up growth and obesity in childhood: prospective 
cohort study. BMJ 2000, 320:967-971. 

32. Coupe B, Grit I, Darmaun D, Parnet P: The timing of "catch-up growth" 
affects metabolism and appetite regulation in male rats born with 
intrauterine growth restriction. Am J Physiol Regul Integr Comp Physiol 
2009, 297(3):R813-R824. 

33. Coupe B, Amarger V, Grit I, Benani A, Parnet P: Nutritional programming 
affects hypothalamic organization and early response to leptin. 
Endocrinology 201 0, 1 51 (2):702-71 3. 

34. Gregor MF, Hotamisligil GS: Inflammatory mechanisms in obesity. Annu Rev 
Immunol 2011,29:415-445. 

35. Li Z, Soloski MJ, Diehl AM: Dietary factors alter hepatic innate immune 
system in mice with nonalcoholic fatty liver disease. Hepatology 2005, 
42:880-885. 

36. Shaikh SR, Mitchell D, Carroll E, Li M, Schneck J, Edidin M: Differential 
effects of a saturated and a monounsaturated fatty acid on MHC class I 
antigen presentation. Scand J Immunol 2008, 68:30-42. 

37. Hua J, Ma X, Webb T, Potter JJ, Oelke M, Li Z: Dietary fatty acids modulate 
antigen presentation to hepatic NKT cells in nonalcoholic fatty liver 
disease. J Lipid Res 2010, 51:1696-1 703. 

38. Zhou B, Wong WH: A bootstrap-based non-parametric ANOVA method 
with applications to factorial microarray data. Stat Sin 201 1, 21:495-514. 



doi:1 0.1 1 86/1 471 -21 64-1 5-49 

Cite this article as: Ellis et al:. Thrifty metabolic programming in rats is 
induced by both maternal undernutrition and postnatal leptin 
treatment, but masked in the presence of both: implications for models 
of developmental programming. BMC Genomics 2014 15:49. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



